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Abstract 



Multivortex dynamics in Manton's Schrodinger-Chern-Simons variant of the Landau- 
Ginzburg model of thin superconductors is studied within a moduli space approximation. 
It is shown that the reduced flow on Mn, the N vortex moduli space, is hamiltonian 
with respect to u>i2, the I? Kahler form on M^v- A purely hamiltonian discussion of 
the conserved momenta associated with the euclidean symmetry of the model is given, 
and it is shown that the euclidean action on (Mjv, 0^2) is not hamiltonian. It is argued 
that the N = 3 flow is integrable in the sense of Liouville. Asymptotic formulae for lo L 2 
and the reduced Hamiltonian for large intervortex separation are conjectured. Using 
these, a qualitative analysis of internal 3- vortex dynamics is given and a spectral stabil- 
ity analysis of certain rotating vortex polygons is performed. Comparison is made with 
the dynamics of classical fluid point vortices and geostrophic vortices. 

1 Introduction 

The Landau-Ginzburg theory of an idealized planar superconductor consists of a complex 
scalar field representing the electron pair condensate, and a £7(1) gauge potential {i = 
1, 2), interacting via the potential energy functional 



Here = <9j0 — iA^ is the gauge covariant derivative, B = d\A 2 — d<iA\ is the magnetic 
field and \i is a coupling constant. The model admits topologically stable, spatially localized 
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solutions (minimals of V) called vortices, the planar analogues of magnetic flux tubes in 3- 
dimensional superconductors. These have finite V, and so have a well-defined winding number 
at infinity, iVeZ, which is interpreted physically as the net vortex number. The case of critical 
coupling, fi — 1, is special because here, given any unordered choice of N complex numbers 
(not necessarily distinct) {z\, . . . , zn}, there exists a minimal of V with iV vortices located 
at z = x\ + ix 2 = z r , r = 1, . . . , N, unique up to gauge equivalence [7j. More precisely the 
Higgs field (f) of the solution has zeroes precisely at the prescribed points z r . Physically, one 
says that static critically coupled vortices exert no net forces on one another. It follows that 
the moduli space of static iV-vortex solutions is = C . Note that the vortex positions 
{z r } are not good global coordinates on Mjy because we must identify solutions which differ 
by permutation of {z r }. Good global coordinates are provided by the coefficients w r of the 
unique monic polynomial in z whose roots are {z r } (counted with multiplicity), that is 

z N + wj/" 1 + --- + w N :={z-z 1 )---{z- z N ). (1.2) 

To introduce dynamics to the theory, one must also define a kinetic energy functional, T, 
and many choices are possible. Manton has advocated the use of a Schrodinger-Chern-Simons 
functional linear in first time derivatives, explicitly, 

T = 7 J (~$ D o<f> - <£D^) + BA o + E 2 A 1 — E\A<2 — A^j d 2 x (1.3) 

where A is the temporal gauge field, Ei = diA — doAi is the electric field and 7 is another 
coupling constant jHj- The resulting dynamics (governed by Lagrangian L = T — V) is first 
order in time and non-dissipative, and is hoped to give a description of vortex dynamics in 
thin superconductors at very low temperatures. The Lagrangian (jl.3|) has also been related 
to effective theories of the fractional quantum Hall effect in ^H] and [23]. Two different 
interpretations for the statistics of the solitons in this model were given in [S] and |15j . 

In order to make progress on the problem of iV-vortex dynamics in this model, one could 
resort to numerical solution of the field equations. Manton [§| has, however, pursued a different 
strategy, namely that of adiabatic approximation. One assumes that at each time the field 
configuration lies in Mjv, the moduli space of static solutions of the \x = 1 system, but that 
the field's position in Mjv may vary slowly with time. This slow variation is assumed to be 
governed by the reduced Lagrangian, that is, T — V restricted to M N . The result is a finite 
dimensional dynamical system, in fact, a hamiltonian flow on M^r, as we shall see, which 
hopefully captures the important features of the full, infinite dimensional field dynamics. Of 
course if \x = 1 exactly, the restricted potential is constant, and there is no dynamics — the 
vortices remain static. On the other hand, if fi differs greatly from 1 then approximating 
multivortex profiles by critical iV-vortices is ill-justified. So the approximation is valid only 
in the near critical regime. This strategy has proved to be extremely successful in the study 
of vortex dynamics in the relativistic version of the Landau-Ginzburg model f7| > and has 
been justified by rigorous error estimates [2*T] . 

The present paper gives an analysis of vortex dynamics within the adiabatic approximation, 
following on from Manton's original paper Pj. Some of the results are anticipated by the 
study of the equivalent model on the two-sphere [TH]. In section |2l for example, we prove 
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that the reduced flow on is hamiltonian with respect to the L 2 Kahler form, a result 
which can be understood from the results on S 2 , at least formally, once the sphere's radius 
is taken to infinity. It is worthwhile giving independent proofs of the relevant results for two 
reasons, however. First, the planar model has far more direct physical significance, both as 
a model of extremely thin superconductors, and of the 3-dimensional case where translation 
symmetry is imposed. The spherical model is useful to study the quantum mechanics of the 
system, but its direct advantage is mathematical convenience (M^v = CP^ is then compact, 
hence the methods of geometric quantization are amenable to more explicit treatment) rather 
than physical applicability. Second, the formulation of the model on S 2 entails considerable 
technical subtleties absent in the planar model. When dealing with the planar model, therefore, 
many of the arguments simplify considerably. 

Neither jH] nor ^3] contains a quantitative discussion of classical vortex dynamics, and 
for good reason — neither the Kahler form nor the Hamiltonian on Mjy are explicitly known. 
In the second half of this paper we will restrict attention to the regime of well-separated 
vortices, where explicit formulae can be inferred, and hence more progress is possible. The 
results turn out to be similar to the dynamics of so-called geostrophic vortices of interest 
in meteorology, and, to a lesser extent, classical fluid point vortices We give a complete 
description of the internal 3-vortex dynamics and analyze the spectral stability properties of 
rigidly rotating regular polygons. 

Lange and Schroers have studied the slow Schrodinger dynamics of ungauged Landau- 
Ginzburg vortices within an adiabatic approximation |8J. For them, there is no concept of 
critical coupling related to a self-duality structure, so no moduli space of static N- vortex solu- 
tions is available. Instead, they define M 2 to be the unstable manifold of the coincident double 
vortex, which they construct via numerical gradient flow. Their emphasis is very different from 
the present work, therefore. They use intensive numerical techniques to study two-vortex dy- 
namics thoroughly, concentrating particularly on the case of overlapping vortices. Indeed, one 
of their main motivations was to provide the first (necessarily numerical) implementation of 
the unstable manifold method for solitons in spatial dimension exceeding one. In contrast, we 
focus on the opposite physical regime (large separation), for N > 3, where heavy numerical 
analysis is not required. 



2 The reduced system 

Let q a , a — 1, ... , 2N, be some choice of local coordinates on Mjv (for example the real and 
imaginary parts of the vortex positions z r G C). Then, since T depends only linearly on time 
derivatives, the Lagrangian restricted to Ma? must take the form 



M, 



A a (q)q a -V(q) (2.1) 



where A a are some functions of the coordinates {q 13 } which may naturally be interpreted as 
the components of a U(l) connexion form A = A a dq a on M^r, since changing A by an exact 
form A i — ► A + d& changes L\m n by an irrelevant total time derivative The Euler-Lagrange 
equations resulting from L\m n are 

dV 

= (2 ' 2) 
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where B = dA is the curvature form of A. Clearly, B is closed. Further, assuming that 
(12. 2|) is a well-defined flow on M^v, B must also be a nondegenerate bilinear form. So we may 
reinterpret B as a symplectic form on M^. Given any function F G C°°(Mjv), let Xp denote 
its symplectic gradient, that is the unique vector field satisfying 

l Xf B = -dF (2.3) 

where i denotes interior multiplication, LxB{Y) := B(X, Y). System ()2.2j) is then hamiltonian 
flow along Xy. 

No explicit formula for B or A exists. However, Manton has shown how A can be related 
to the analytic properties of the Higgs field <ft near its zeroes. Let us assume that the vortex 
positions z r are all distinct, and use these as local complex coordinates on M^r- Then it is 
known that in a neighbourhood of z r 

log |0(2)| 2 = log \z - z r \ 2 + a r + -& r (z - z r ) + -b r (z — z r ) H (2.4) 

where a r and b r are some unknown functions of the vortex positions z s , a r real, b r complex. 
We have used here the conventions of [llj. Samols [17J, adapting earlier work of Strachan 
[2*U] . showed that the L 2 metric on M^v could be written entirely in terms of the unknown 
coefficients b r , namely 

N / F)h\ 

(2.5) 



( db s \ 

r,s=l r / 



This is a hermitian metric on TMjv which can be shown to be Kahler ^7j. Using similar 
techniques, Manton has obtained jU] a similar formula for the connexion form A; in our 
notation, 

N r - 1_ 1 

{b r + -z r )dz r — (b r + -z r )dz r 



A = ITX^j 



r=l 



(2.6) 



We will now establish the following proposition: 



Proposition 1 The connexion A has curvature form B = —l^ioip. where ujl? is the Kahler 
form corresponding to g^ . 

Proof: We will prove the formula on M^An, where A^v is the measure zero set on which 
vortices coincide, and appeal to smoothness. Let us define the (0, l)-form 



b:=^2b r dz r . (2.7) 



Hermiticity of gip. and (|2.5jl imply 

db. db 



TT— = 7T- =► 8b = -db. (2.1 

OZ r OZc 
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Note that the connexion A and lOip. may be written, using ()2.7j) and (|2.8|) 



A 



27T7 



- S ^^(z r dz r — z r dz r ) + b — b 



U L 2 



ITi 

Y 



yj dz r A dz r + db — db 



(2.9) 



It follows that 



2uo L 2 + -dA = in(db -db + db-db) = m{db -db). 
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(2.10) 



Now ddb = —ddb = d 2 b = by ()2.8|) . so db is an antiholomorphic (0, 2)-form. Hence its 
component functions are antiholomorphic functions on Mtv\Atv- But it is known J7| that the 
coefficients b r , and their derivatives, decay exponentially fast at large vortex separation. Hence 
these antiholomorphic functions must in fact vanish. So db = = db, and the proposition 
follows. □ 
For future reference, we note that db = implies that 



db r 

dz R 



dK 

dz r 



(2.11) 



for all r, s. 

So the adiabatic approximation to iV-vortex dynamics is hamiltonian flow on (M^jUJ^ 
with Hamiltonian 



H 



1 



-V\ 



Mr 



(2.12) 



Henceforth Xp will denote the symplectic gradient of F with respect to ujl2, rather than B. 



3 Conservation laws 

We may now give a hamiltonian discussion of the conservation laws discovered by Manton 
and Nasir jlOj . The natural action of the euclidean group E(2) = U(l) x C on the physical 
plane C induces a S(2)-action on M N by (e t9 ,c) : {z r } \— > {e te z r + c} (strictly speaking, this 
only defines the action on Mtv\Ajv, but one can use (jl.2|) to deduce a well defined action on 
the global coordinates w r ). This action is manifestly holomorphic. It is also isometric, since 
it leaves L 2 -norms invariant. Hence the £'(2)-action is symplectic. We would like to construct 
a moment map fi : M^v — > e(2)* for this action and identify its components with respect to a 
natural basis for e(2)* as the conserved momenta of the system, just as it has been done for 
the symplectic SO (3)-action in the S 2 model [T5] . Unfortunately, no such moment map exists. 

Certainly, given any X G e(2) = T {lfl) E(2), the induced vector field X s G T(TM N ) is 
hamiltonian, because all symplectic vector fields on are hamiltonian (Y symplectic implies 
lyojip. is closed; it is also exact since if^M^v) = 0, whence Y is the symplectic gradient of some 
smooth function). Using the terminology of ^2], the i?(2)-action on is weakly hamiltonian. 



5 



So functions Pj : Mjv — ► R , j = 0, 1, 2 satisfying — <iPj = ix^t?- must exist, where 

9" ( d d \ 



<9ci ^— ' \cte r 9l 

generate rotations and translations in M^r. In fact, these functions are unique up to addi- 
tive constants. The problem is that the Pj cannot be assembled into the components of an 
equivariant map Mjv — > c(2)*. In other words, we have: 

Proposition 2 The symplectic E(2)-action on is not hamiltonian. 

Proof: Assume, to the contrary, that a moment map \i exists. Each M^v contains a symplectic 
submanifold (C, 1I ^-dZ A dZ), namely the P(2)-orbit of any configuration of N coincident 
vortices (here Z denotes the common vortex position). The restriction (i of /x to this orbit 
defines a moment map for the standard action of E(2) on C. Since fx is equivariant, and the 
5*0(2) action fixes 0, 

fi(0) = K(e l6 ,0) -0) = Ad* (e , efl) m (3.4) 

where Ad* denotes the coadjoint action of E(2) on e(2)*. Hence fi(0) = ad9 for some real 
constant a. But this completely determines fi by transitivity of the translation action: 

ft(Z) = A((l, Z)-0) = Ad\ l z) ade = ad6. (3.5) 

Thus, any equivariant map C — > c(2)* is constant, and cannot generate a nontrivial group 
action. □ 
Notwithstanding the lack of a moment map, generating functions Pj for the Xj do exist, 
and one expects, given Manton and Nasir's Lagrangian analysis ^H], they may be written 
locally (on Mjv\Ajv) as 

P = - ^(|<2v| 2 + b r z r + b r z r ) (3.6) 

r 

Pi = *~I>r-*r) (3-7) 

r 

P2 = l/Z^+^r)- (3-8) 

r 

Note that P\^P% are easily globalized since Pi + iP<i = —mvj\ in the global coordinates. A 
smooth globalization of Pq must exist, but no obvious formula suggests itself. 

It is a routine exercise to verify that dPj(Y) = u^iY, Xj) for all Y G T(TMn) as required, 
if one makes use of the identities 

J> = (3.9) 

r 

^2{b r z r -b r z r ) = 0. (3.10) 
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Identity (|3.9|) is proved in J7|, while (|3.10|) is proved for the spherical model in [15J. For the 
sake of completeness, we will derive it directly for the planar model. 

Let (p G Mat be the Higgs field of the static solution with zeroes at {z r }, and let 4> e <E Mn 
be the field obtained from this by rotating all the vortex positions Then clearly 

4> (e l9 z) = 4>(z), and hence ()2.4j) implies that 

b r ({e id z r })=e id b r ({z r }) (3.11) 

for all r. It follows that 

b r = ib r , and X [6 r ] = — ib r , (3-12) 
since X is real. But using our local expression for X , (|3.1|h we see that 

where we have used (|2.8j) and (|2.11|) . Comparing (|3.13|) and its complex conjugate with (|3.12|) 
one sees that ^ r {b r ~z r —b r z r ) is constant. That the constant vanishes follows from consideration 
of the limit where all N vortices coincide. In fact, constancy of the sum, rather than vanishing, 
is sufficient to prove that Xp = X . 

Our Hamiltonian H = — t^^Imat is manifestly invariant under the _E(2)-action, so 

{P^H} = uo^X^Xh) = dH(X 3 ) = Xj[H] = (3.14) 

and the momenta Pj are conserved by the flow of H . 

The Lie subalgebra generated by the Pj in the Poisson algebra C°°(M7v) is determined by 
the relations 

{Pj,P } = e jk P k , j = 1,2 (3.15) 
{P 1 ,P 2 } = Ntc, (3.16) 

where denotes the antisymmetric tensor with ei2 = +1. So we can see that it is a central 
extension e(2) © R of the euclidean algebra. The nonvanishing of the Poisson bracket (J3.16|) 
is consistent with a field theory calculation of Hassaine et al. [6]. In our context, given 
that [Xi,X2] = 0, it provides another proof of Proposition^ In fact, it implies that the 
_E(2)-action on Mjv corresponds to a nonzero 2-cocycle in the Lie algebra cohomology group 
iJ 2 (e(2);M) = R; this group parametrizes obstructions of weakly hamiltonian £'(2)-actions to 
be hamiltonian [T2] . 

It is worthwhile to look at the content of Proposition |21 in the context of the results on the 
spherical model jT^j. In the latter, spatial isometries are described by the group SO (3), which 
is simple. Therefore, if 2 (so(3);R) = and the S'0(3)-action on Mjv = CP^ is necessarily 
hamiltonian. Dualizing the corresponding moment map, we obtain an isomorphism between 
so (3) and the Lie algebra of conserved angular momenta. The large radius limit R — >• oo that 
relates the spherical model to the planar model determines a contraction of Lie algebras 

so(3) -> e(2) (3.17) 



MK] = -to 
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with parameter ~. This is the well-known contraction of so (3) to the graded Lie algebra 
associated to the filtration jlj 

C so(2) C 50(3), (3.18) 

where the choice of the so (2) subalgebra corresponds to fixing an axis through the centre 
of the sphere. We can realize this contraction as a deformation of Lie algebras of vector 
fields, using the isometry actions on each M^. At the level of the Poisson algebras, however, 
no homomorphism of Lie algebras from e(2) to C°°(M7v) is available by Proposition so 
it is not surprising that the closure of the generators Pj does not yield e(2) in the planar 
model. A completely analogous discussion can be given for the zero-curvature limit relating 
the hyperbolic plane to the euclidean plane. 

To end this section, we would like to emphasize that the fundamental reason underlying 
the nonhamiltonian action of the isometry group is the geometry of the domain of the model, 
rather than the dynamical nature of the model itself. Indeed, the proof of Proposition El relies 
crucially on the fact that the standard action of E(2) on the plane is not hamiltonian, so 
one would expect this to be a common feature of symplectic moduli spaces of planar solitons. 
By contrast, the equivalent symmetry groups for the sphere and hyperbolic plane, SO (3) 
and 5*0(2, 1) respectively, do have hamiltonian actions on the corresponding vortex moduli 
spaces. The lack of a moment map means that the standard Marsden-Weinstein technique 
for reducing the dynamics to low-dimensional symplectic quotients is not available to us. The 
Poisson noncommutativity of the components of the linear momentum (|3.1fij) is the "classical 
anomaly" which provides the obstruction. It is an expression of the fact that two of the real 
coordinates on phase space (essentially, the real and imaginary parts of Wi) are conjugate for 
the relevant symplectic structure, and should not be interpreted as a breaking of translational 
invariance. 



4 Well- separated vortices 



Although the moduli space approximation has simplified the dynamical iV-vortex problem 
in principle, we are still faced with the problem that neither the Kahler form loui nor the 
Hamiltonian H = — ^-V are explicitly known. There is one physical regime in which explicit 
progress is possible, namely when all vortices are well separated from one another. In this 
situation, an asymptotic formula for the coefficients b r , and hence uo L 2 is available jTT]. By 
translation invariance, the formulae involve only Zy a • Zy Zg • They are 



. 7T 



0>£2 = l- 



b r (zi, ...,z N ) 



2tt 



dz r A dz r — ^— ^ K (\z rs \)dz rs A dz r 

r r s^r 



(4.1) 
(4.2) 



where K n denotes the modified Bessel's function of the second kind, and q is an unknown real 
constant which may be interpreted as the scalar monopole charge or magnetic dipole moment 
of a single vortex at critical coupling. Based on a string theoretic duality argument, Tong 
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has conjectured that q = — 27r84, and this value is consistent with numerics [221 The 
formulae are believed to be correct up to linear order in exponentially small quantities (note 
K n (p) is exponentially small for large p). They have not been proved rigorously, but they 
have been derived via two independent physical arguments, one using matched asymptotic 
expansions and the other appealing to a point source formalism, and there is a good fit to 
known numerics. Similarly derived formulae for monopole moduli spaces have subsequently 
been proved rigorously, so one has good grounds for confidence. 

It is not so straightforward to obtain an asymptotic formula for H. It is known [§] that 

V\ Mn = Nn + ^± J {I- |0| 2 ) 2 d 2 x, (4.3) 

but no obvious method of estimating the integral for well-separated vortices suggests itself. 
We shall instead proceed by an indirect physical argument. For p close to 1, one expects that 
V\u N is close to V£, the Landau-Ginzburg energy of iV well-separated coupling p vortices 
held at positions {z r }, and for this an asymptotic formula has been found |19j : 

VS^EEwCk.!) (4.4) 

r sj^r 

where U is the asymptotic 2-vortex interaction potential, 

U{p) := Um{pfK Q {p) - q{p) 2 K {pp)l (4.5) 
Ztt 

and we have discarded an irrelevant constant (N times the vortex mass). The factor of \ is 
included in (|4.4j) since the sum is over all ordered distinct vortex pairs. Note that the scalar 
monopole charge q(p) and magnetic dipole moment m(p) vary with the coupling, and are not 
equal except when p — 1. Note also that, at p — 1, 11 = 0, as one expects. Now let p = l + 5p, 
Sp small, and assume that, at sufficiently large vortex separation, V\m n = V£ to leading order 
in Sp. Again discarding additive constants, one is led to conjecture that 

r sj^r 

U{p) := f-(p* - 1) {p Kl (p) - vK (p)} , (4.6) 

where v = 2(g'(l) — m'{l))/q w 2.7060, using the tabulated data for m(p) and q(p) in [T9j . 
Note that equation ()4.3)1 is an exact equality which is not restricted to the near critical (p « 1) 
regime: the Landau-Ginzburg energy restricted to M^v is precisely proportional to p 2 — 1. Our 
asymptotic formula for it (in the well-separated regime) is obtained by assuming approximate 
equality with V£, an assumption which can only hold for p close to 1; but the formula ()4.6|) 
itself is not restricted to the near-critical regime. Of course, the validity of the whole moduli 
space approximation becomes questionable for \p — 1| large, but this is a separate issue. 

One way to test our conjecture is to deduce from it an asymptotic formula for the integral 

/shah(p) = |(l-|0r) 2 rf 2 x (4.7) 
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2 4 6 8 10 12 14 

P 

Figure 1: The quantity J*(l — |0| 2 ) 2 <i 2 x for a vortex pair, as a function of the vortex separation 
p. Solid line and crosses: Shah's numerical data, taken from JH]/ dashed line: our conjectured 
asymptotic formula, as in equation (14. 8|) . 

for the case of critically coupled N = 2 static vortices with intervortex separation p (note that 
E(2) symmetry implies the integral may only depend on p in this case). The function /shah 
has been evaluated numerically by Shah JHI- Formula (|4.6j) implies that 

/shah(p) - /shah(oo) = ^-[pK x {p) - vK Q {p)\ (4.8) 

7T 

asymptotically at large p. A comparison of Shah's numerical data with our conjectured asymp- 
totic form is given in figure The fit for p > 4 is very good. 

The iV-vortex equations of motion in the asymptotic regime are obtained by extracting the 
leading asymptotic term in the symplectic gradient of H. Since H is already exponentially 
small, an exponential correction to uj^s makes no contribution at this order, so we can actually 
set U£2 — nu , where 

uj = % - ^ dz r A dz r (4.9) 

r 

is the canonical symplectic form on C^. We find, to leading order, 

z r = iY j F{\z rs \)z rs , F(p) := (4.10) 

s^r 

Note that system ()4.10j) is similar to the equations of motion for a system of identical fluid 
point vortices, or geostrophic vortices, of vorticity I/7, which would result from replacing U 
by 

Wfluid(p) = logp, W g eo(p) = K (p) (4.11) 
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respectively [HUB]- The system with £4 U id is particularly well studied. We shall see that our 
system behaves rather more like the geostrophic vortex system, which has not been so heavily 
studied, though there are still significant differences. 

Two-vortex dynamics in the moduli space approximation is almost trivial: the vortices 
orbit their centre of mass at constant speed and separation. This is a special case of the 
rotating polygon solutions whose stability properties we will analyze in section We turn 
now to a discussion of 3- vortex dynamics. 



5 The dynamics of three vortices 

Novikov has given a thorough treatment of the internal dynamics of identical vortex triples 
interacting via Wfl ui d in [T3] . His method can be adapted readily to deal with our system in its 
asymptotic form, (|4.10jl . The basic idea is to identify trajectories in the internal phase space 
of the system with level curves of the quantities conserved by flow (|4.1()j) . Since Novikov's 
method relies on exploiting the conservation laws enjoyed by the system, it is more satisfactory 
to apply it directly to the momenta conserved by the full flow, in asymptotic form, rather than 
the momenta conserved by the asymptotic flow ([4.10)1 . and this is how we shall proceed. Of 
course, either approach yields the same results. 

To begin with, let N be general. Let us define the centroid of an iV-vortex configuration 

by 

r 

and note that 



2JV 



Since Z and Pq are conserved, it follows that 

27V 

7T 



^ ( N(b r z r + b r z r ) - ^ ^2 \ z n 

r \ sj^r 



(5.3) 



is conserved also. The corresponding hamiltonian vector field can be readily computed as 

X Q = ™X --^(2P 1 X 1 + 2P 2 X 2 ) 

and this describes rigid rotations of vortex configurations about their centroids. Using (|2.8J) . 
(gHH) and (prmjl . one finds 

{P J ,Q} = 0, j = 1,2. (5.5) 
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For a moment, let us consider N = 3. Equation ()5.5|) implies that the set of three con- 
served quantities {Pi, Q, H} is in involution. If the corresponding hamiltonian vector fields are 
linearly independent at all points of an open dense subset M3 \ S C M3, the dynamics of three 
vortices is Liouville integrable jH] on the six-dimensional phase space M3. This manifold is then 
foliated by invariant tori some of them possibly degenerate. Clearly, X\ is independent 
of both Xq and X# everywhere in M^. Hence, the question of Liouville integrability reduces 
to the condition of linear independence of Xq and X#. Our intuition about the dynamics of 
N = 2 vortices [§| leads us to claim that this condition is generic on M 3 . Actually, we expect 
Xq and X# to be linearly dependent exactly at the points of a subset S C M 3 consisting of 
configurations with special symmetry: z r either at the vertices of an equilateral triangle, at 
the ends and midpoint of a line segment, or on the locus of coincidence A 3 ; notice that all 
these conditions are algebraic in the coordinates z r (and their complex conjugates), so that 
S is indeed closed with dense complement. We shall see in the following that this picture is 
consistent with what we can learn about the dynamics of 3-vortices in the asymptotic regime 
of large separation. The configurations in S that correspond to well-separated 3-vortices are 
special cases of the vortex polygons that we shall discuss in section 

When all \z rs \ are large, the angular momentum becomes 



Pn 



7T 



r r s^r 



whereas 



«4EE 



Nq 2 

\Kx{\ 



7T" 



(5.6) 



(5.7) 



Note that, unlike Z and P , Q and H depend only on the relative separations of the vortices 
\z rs \, to this order. 

If we now set iV = 3 and define x = (| z 2 z |, \z$\ |, | -212!), we see that x(t) is confined to a 
level curve of (H, Q) in R 3 , where 



H 

Q 



U(xt) +U{x 2 ) +U(x 3 ), 



Nq 2 



g(xx) + g(x 2 ) +g(x 3 ), g(x) := x H —xK^x) 



7T 



(5.8) 



Our strategy is, then, to construct such level curves. This almost completely determines 
the internal dynamics of the triple (meaning the dynamics up to rigid rotations), up to 
reparametrization of time. The only thing not determined is a discrete variable, the ori- 
entation of the triple. Two orientations are possible: traversing the triangle with vertices at 
z r clockwise can order the vertices 123 or 132. 

To construct a level curve of (H, Q), we may choose some initial data x(0) = xo and then 
solve the first order system 

x(s) = (5.9) 
\VH x VQ| 

to yield a level curve in arc-length parametrization (|x(s)| = 1). Since we work only to linear 
order in exponentially small quantities, and H is small, we may keep only the leading term in 
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y =0 

J 3 

Figure 2: Level curves of H Q for Q = 441. Here Hi = H® m and H 2 = H^a- The dashed 
curves represent the colinearity boundaries, which intersect at coincidence points. 

Q, that is, take 

<? = |x| 2 (5.10) 

in (|5.9|) . So level sets of Q may be approximated by concentric spheres of radius in M 3 , 
and the level curves of (H, Q) are simply the intersection curves between level sets of H and 
these spheres. Alternatively, given a radius Q^, we may seek level curves of H Q : S 2 — > R, 

^(y) = W(g^x) + W(Q5y 2 ) + W(Q*2/ 3 ), (5.11) 

where x = Q^y and 5* 2 is the unit 2-sphere in M 3 . It is straightforward to generate contour 
plots of H® for given Q within any convenient coordinate system on S 2 , and hence deduce 
the corresponding internal trajectories. Note that Q determines the root-mean-square of the 
triangle side lengths, xi,X2,xs. It turns out that the level curves are not very sensitive to Q, 
provided we choose it to be fairly large (which it should be, in the asymptotic regime), so let 
us imagine that Q has been fixed. 

Not every y G S 2 represents a valid vortex triple. Each must be non-negative, and 
since they represent the side lengths of a euclidean triangle, none may exceed the sum of 
the other two. So y(t) is confined to a triangular region in S 2 , bounded by great-circular 
arcs (the intersection of the sphere with the colinearity planes y\ = y 2 + 2/3, y 2 = 2/3 + 2/1, 
2/3 = 2/i + 2/2), and centred on y = (1, 1, l)/3a. This central point represents an equilateral 
triple. The orthogonal projection of this region onto the tangent plane T yo S 2 is depicted 
in figure 121 The boundaries project to congruent elliptic arcs. Any point on the boundary 
represents a colinear vortex triple, and the boundaries intersect at coincidence points (where 
a vortex pair coalesces). Since our analysis is restricted to the well-separated regime, level 
curves which approach these corners closely are of suspect validity. 
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Figure |21 depicts level curves for the choice Q = 441 and several values of H. Our discussion 
is qualitatively similar to Novikov's [Tl], so will be kept brief. The minimum possible value of 
H® is 

^ m = ^(yo) = 3W( V / g/3) (5.12) 

attained by the equilateral triple. The level "curve" for this value consists of the single point 
yo, so the internal configuration of this triple remains constant. In fact the triangle rotates 
about its centre at constant speed, as will be discussed in section |3 Increasing H Q , we 
obtain closed curves centred on y . In this low energy regime, the shape of the triangle varies 
periodically in time as the triangle spins, but the triangle's orientation remains constant. The 
periods of shape change and rotation are generically incommensurate. Clearly, these solutions 
are all stable. 

Increasing H® further, the level curve eventually becomes tangent to the colinearity bound- 
aries, at the critical value 

H2 it = l,2)/6*) = 2W(y^/6) +W(2y / g/6). (5.13) 

This level curve should not be interpreted as a closed orbit. Rather, we have three fixed 
points, (1, l,2)/\/6, (1,2, l)/y/E, (2, 1, l)/y/E, joined by 3 heteroclinic orbits. The fixed points 
represent colinear triples with the middle vortex equidistant from the outer vortices, a special 
case of the filled rotating polygons which will be discussed in section |3 This is stationary (but 
not static) and clearly unstable: small perturbations lead to divergent dynamical regimes. 

For > H2 it , the level curve splits into 3 disjoint pieces. The actual trajectory stays on 
one segment of this curve - it reverses direction each time it hits a colinearity boundary, with a 
corresponding change in the triangle's orientation. At these points, the triangle collapses and 
turns itself inside out, so for larger energies, the triple's internal motion consists of a periodic 
"flip-flopping". As the energy grows unbounded, the trajectory approaches a coincidence 
point, and hence escapes the presumed region of validity of our approximation. 

It is interesting to note that every trajectory intersects the isosceles lines (the straight 
lines containing y and the coincidence points), so every 3- vortex motion passes, at least once, 
through an isosceles configuration. 

6 Rigidly rotating polygons 

If N vortices are placed at the vertices of a regular iV-gon, they will rotate at constant 
speed about their centroid, which remains fixed This behaviour is familiar from both the 
fluid and geostrophic vortex systems. Given the asymptotic equations of motion (|4.10p . it is 
straightforward to derive frequency-radius relations and spectral stability properties of these 
solutions, valid for large polygon side length. Note that a change in the coupling constants /i, 

7 may always be absorbed into a rescaling and/or reversal of time, so the stability properties 
will be independent of our choice of parameters. 

It is convenient to transform to a co-rotating frame by defining new coordinates y r (t) G C, 

z r (t) =: y r (t)e mt (6.1) 
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for some fixed choice of frequency f2 G R, so that frequency solutions are static in the new 
frame. System (j4.10|) in these coordinates is 

y r = i^2F(\y TS \)y rs -iSly r , (6.2) 

s^r 

where y rs := y r — y s . Note that ()6.2|) is hamiltonian flow on (C N ,u> ) with respect to the 
modified Hamiltonian 

Let A := e 2m l N . Then y r = A r cr, a G (0, oo), is a static solution of ()6.2j) provided 

N-l 

ft = ^F(|l-A>). (6.4) 
i=i 

This is the frequency-radius relation for a rotating iV-gon. 

We turn now to stability properties. Let y G be a fixed point of the co-rotating flow. 
Then we may define the linearization of the flow about y, 

A y : T y C N -> T y C N , A y : X ^ V x X Hn , (6.5) 

where V is any afline connexion on C N , and X# n is the symplectic gradient of Hq. The fixed 
point y is spectrally stable if A y has no eigenvalues with positive real part. Now, since X/f^ is 
symplectic, it follows that A y G aut(T y C N , u ), that is, 

u {A y X, Y) + u {X, A y Y) = (6.6) 

for all X, Y G T y C N . It follows that, if r/ is an eigenvalue of A^, so are —rj, r\ and — r\. In order 
that y be spectrally stable, therefore, specA^ must be purely imaginary. 

In the case of the rotating iV-gon, y = y = (Act, AV, . . . , ct), one can find spec A^ explicitly 
due to the cyclic symmetry of the configuration jT3|. Let tt : Z — > Zjv be the group homomor- 
phism with 7r(l) = 1. Then a N x N real matrix A is right circulant if there exists A G M. N 
such that 

Aij = A^-j+i) (6.7) 

for all i, j. In order to exploit this symmetry, it is convenient to define a twisted complex basis 
e(j) G C^, i = 1, . . . , N, such that 

e (i)j = A% (6.8) 

there being no summation implied. If we write down the matrix representing Ay with respect 
to the associated basis for M. N , that is, {Ree(i), . . . , Ree( n ), Ime(i), . . . , Ime( n )}, we find that 



Ay 



A B 
C A 



(6.9) 
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where each of the real N x N blocks is right circulant. The explicit formulae for A, B, C G M. N 
are rather complicated, and are presented in the appendix. Of more interest are the following 
symmetries possessed by these coefficients: 



An+2-s = 


-A s 


Bn+2-s = 


B s 


Cn+2-s = 


C s 



. , N. (6.10) 



Right circulancy is a powerful symmetry because the eigenvectors of a right circulant matrix 
are independent of its entries. For each j G {1, . . . , iV}, let 



x 



%\ j ,\ 2j ,...,\ (N - 1)j ) G (6.11) 



Then x 3 - is simultaneously an eigenvector of A, B and C, with eigenvalues 

a j = A-x j , (3 j = B-x j , 7j = C'-x i . (6.12) 

From (|6.10j) . (jfi.llj) and ()6.12|) . one sees that aj is purely imaginary, while (3j and jj are real. 
Note that Xj = xjy^j for all 1 < j < N — 1, and that Xjv = x^. It follows that 

a N -j = aj = -aj, (3 N -j = = (3j, ^ N -j = %• = lj, (6.13) 

for 1 < j < N — 1, and that «at = = —&n, so = 0. Short calculations using the 
detailed forms of B, C (see appendix) establish that (3n = and 71 = —f3i. 

Now, relative to the basis {x^ , xj", . . . , x^-, x^}, where x^ := (xj, ±Xj) G C 27V , is block 
diagonal, with 2x2 blocks 



A, 



ay + i(ft + 7i) l^i + Ti) 



(6.14) 



It follows that the eigenvalues of Ay are 



'// = <*i ± V 7 ^- (6-15) 
Owing to the symmetry properties of aj, (3j, 7^, (jEHHJ), we see that 

= -vf, (6.16) 

so that eigenvalues generically come in quartets. Note also that 7/^ = (since a at = (3n = 0) 
and that rji = —i]n-i = These eigenvalues originate from the E{2) symmetry enjoyed by 
system f)6.2jl . as will be shown shortly. 

Since aj G iM. for all j, y is spectrally stable provided that fij^j < for all j. But (3\ = —71, 
and Pw-jlN-j = Pjlj for 2 < j < N — 1, so there are in fact only k — 1 stability criteria, 
where N — 2k or N — 2k + 1, namely 

/Vtt<0, J = 2,...,fc. (6.17) 
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N 


2 3 4 5 5 7 >8 


gauged 
fluid 
geostrophic 


11110 
111111 
11110 



Table 1: Comparison of spectral stability properties for rotating vortex N-gons of various types, 
=unstable, 1 =stable. In the gauged and geostrophic cases, the entries refer to stability at 
large polygon radius a. There may be windows of anomalous (in)stability at small a, but these 
lie outside the regime of validity of the calculations. Entries for geostrophic and fluid vortices 
are taken from [T3] . 

For iV = 2, 3 stability is automatic, for N = 4, 5 there is one criterion, for N = 6,7 two 
criteria, and so on. 

In the first nontrivial case, N = 4, A = i and the single stability criterion is 

(/3272)0) = -8V2a 2 F'(V2a)F'(2a) < 0. (6.18) 

Since F has no critical points in the well-separated (large a) regime, the rotating square is 
spectrally stable for all a sufficiently large. As N increases, the expressions for f3jjj become 
increasingly complicated so that it is not feasible to check their sign by hand. It is straight- 
forward to check the criteria graphically by plotting the sign of (3j^j against a , however. The 
results are summarized and compared with the previously studied cases of fluid and geostrophic 
vortices in tabled 

We may give a similar analysis of the case of A^ — 1 vortices on the vertices of a regular 
polygon orbiting another vortex at the polygon's centre. Let A = e 27r4 /( Af - 1 ) # Then y = 
(0, Act, A 2 cr, . . . , cr) G is a fixed point of the flow (jfi.2j) . provided that 

N-2 

Q = F(a) +J2 F (\ 1 ~ - xj )- ( 6 - 19 ) 

j'=i 

Once again the linearization G out(T^C 7V , uj ), so y is spectrally stable if and only if 
specA^ C iM.. The presence of the centre vortex destroys the right circulant symmetry of 
A, however, so we must resort to a numerical algorithm to generate specA^. By plotting 
max{|Re^| : rj G specA^} against a, we can determine which of the (N — l)-gons are spec- 
trally stable in the well separated regime. The results, which were generated using Matlab's 
eigenvalue finder, are summarized in table El 

These results have been derived in the approximation of large a. In the full adiabatic 
approximation, they will receive corrections due to the subleading terms in lou^ and H, and 
we are assuming that these corrections will be very small when a is large. What reason do 
we have to believe that the seemingly delicate property of spectral stability, spec A C iM, will 
not be destroyed by these perturbations? The spectra found exactly (for A^), or numerically 
(for Ay), had no accidental degeneracies for a sufficiently large, in the stable cases. All the 
eigenvalues were simple, with the exception of rj — 0, which we will shortly return to. No 
simple imaginary eigenvalue can be perturbed off the imaginary axis because, due to the 
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N- 1 


23455789 > 10 


gauged 
fluid 
geostrophic 


111110 
1111111 
1 1 1 1 



Table 2: Comparison of spectral stability properties for vortex (N — l)-gons rotating about a 
single central vortex for various vortex types, key and comments as for table 

reflexion symmetries enjoyed by spec A, it would have to split in two, which is impossible by 
conservation of multiplicity. Further, the eigenvalue, which turns out to have multiplicity 2, 
is fixed at by symmetry considerations, as we shall now see. 

We are free to transform to the co-rotating frame in the full adiabatic flow, by redefining 
our Hamiltonian 

H i-> H n = H- ttP . (6.20) 

Then y and y will still be fixed points of X Hn , provided a is chosen suitably. In each case (y — y 
or y) we may study the spectrum of the linearization A y : X i— > Vx^h q G out^Mjv-, ^i 2 ), 
which will be a slightly perturbed version of the asymptotic spectrum already considered. 
Recall that X , Xi, X 2 denote the vector fields generating rotations and translations (|3.1|) - 
(13. 3|) . Now for any fixed point y <E M N , 

\ y X 3 = V x X Hn = [X S) X Ha ] = -n[X S) X ] (6.21) 

since [Xj,X H ] = X{p #} — by the E{2) symmetry. Hence 

K y x Q = o, A y X! = -nx 2 , A y x 2 = nx u (6.22) 

whence it follows that {0,iQ, —ifl} C specA y . Hence, the double eigenvalue is fixed under 
any perturbation maintaining E(2) symmetry, so we conclude that the simple stability analysis 
given above is structurally stable at large a. 

7 Conclusion 

In this paper we have confirmed the expectation raised by ^H] that the adiabatic approximation 
to A^-vortex dynamics in Manton's model of a planar superconductor amounts to a natural 
hamiltonian flow on (Mjv,^l 2 )- A hamiltonian account of the conserved momenta descending 
from the symplectic but nonhamiltonian action of E(2) on Mjy has been given, and based on it 
we argued that the dynamics of three vortices is integrable in the sense of Liouville. We have 
derived an asymptotic formula for the 2-vortex interaction potential close to critical coupling, 
namely 

U{p) = |V - l)[pKiO) - uK (j>)], q = -2tt83, v « 2.7060, (7.1) 

and used this to analyze internal 3-vortex dynamics at large separation. We have studied 
the spectral stability properties of rigidly rotating vortex polygons, both with and without a 
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central vortex, and compared our results with the more familiar cases of fluid and geostrophic 
vortices. 

Now that we have some quantitative dynamical predictions from the adiabatic approxi- 
mation, it would be interesting to test its validity numerically. There are really two separate 
issues. First, does the hamiltonian flow give a good account of low-energy vortex dynamics 
near critical coupling? Second, do our asymptotic formulae for U and u L 2 give a good ap- 
proximation to this flow at large separation? Both questions may be addressed by numerical 
analysis of 2-vortex dynamics. According to the adiabatic approximation, two vortices orbit 
one another at constant radius and angular frequency indefinitely. Is this behaviour seen in 
numerical simulations of the full field theory? One would expect the orbiting pair to gently 
radiate energy and hence slowly drift either together or apart (depending on the choice of p, 
7), but hopefully this happens on a much longer time scale than the period of their orbit. If 
their separation p is large, then their angular frequency should be close to 



This is very different from the inverse square law found for ungauged Landau-Ginzburg vor- 
tices 0. It would be interesting to test this formula numerically in terms of both its p and p 
dependence. 
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Appendix 

The simplest way to construct the matrix representing Ay with respect to the twisted basis eu) 
defined in equation ()6.8|) is to incorporate the twisting into the co-rotating coordinate system, 
by defining x r (t) such that y r (t) = \ r x r (t). The equations of motion (|6.2|) become 



Then x — (a, a, . . . , a) is a fixed point of (jAlj) provided (J6.4j) holds. The linearization of (jAlj) 
about x is 



5x r = iR(a)5x t - iS(a)5x r - i ^ A s - r Q(|l - \ s - r \a)5x s + i J^Pfll ~ A s -»5x s , (A2) 




(7.2) 




(Al) 



where 




F(a)+P(a), 
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JV— 1 



R(a) := ^F(|l-A>) + P(|l-A>)], 
3=1 

JV-1 

S(a) := ^A'P(|1-A>). 



(A3) 



3=1 



Note that P,Q,R,S are all real. Decomposing (jA2|) into real and imaginary parts, with 
6x r = 6x1 + ibxi, one finds that 



6x\ 



6X]y 

6xj 



6x 2 N 



A B 
C A 



l n 



6x\ 



8x 1 N 
6x\ 



6x 2 N j 



(A4) 



where A, P, C are the N x N right circulant matrices generated by A, P, C G M. N : 

s = 2, 



A 1 = A^Im^'-^Qfll- A s -V) 

Bi = -P(a) - 5(a) 4 = Re (A S - X )Q(|1 - A^ja) + P(|l - A S ~V) 



s = 2, 



d = iZ(o-) - S(<r) C s = -Re(A s - 1 )g(|l-A s -V) + P(|l-A s - 1 | ( x) s = 2, 



N 
N 
N 



(A5) 

The matrix in equation (|A4|) is A^ in this basis. The symmetries of the components of 
A,B,C claimed in (j6.1U|) follow immediately, once we note that A 
that = X s ' 1 w hen 2 < s < N. Note that 



JV 



1 and A" 1 = A 



so 



JV 



p. 



N 



J^B.s = -R(a) - S{a) + [ Re (A S_1 )Q(|1 - A S "V) + P(\l ~ ^W)] 



8=1 

N-l 



s=2 



as previously claimed. Note also that 

JV-l 

A + 7i = Pi + 4 + ^2(B j+ i + C j+l )X j 

3=1 

N-l 

= -25( ( j) + 2^A- ? P(|1-A>) = 
so that 71 = —fli, as was claimed. 



(A6) 



3=1 



(A7) 
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